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Abstract—A simple and efficient procedure is 
proposed to predict the performance of large finite 
arrays from the measurements of a small array of 
similar lattice. The procedure depends on_ the 
knowledge of the mutual admittance matrix of the 
small array from which the mutual admittance 
matrix of the larger array is constructed by 
ignoring the mutual coupling between the elements that 
are external to the domain of the small array. Once 
the admittance matrix has been constructed for the 
large array, active input impedances of the array 
elements under different scenarios can be predicted 
easily. Also, from the active voltage terminal of the 
elements, an array factor including the mutual 
coupling can be computed to predict the radiation 
patterns as well as the array gain. Here, an air 
microstrip patch antenna excited by a hook-shaped probe 
is considered as an array element. The results obtained 
using the present method is verified by the full-wave 
numerical analysis of the large array. The results 
obtained serve to verify the concept. 


Index Terms—Finite Array, Large Array, Microstrip 
Antenna, Phased Array. 


I. INTRODUCTION 


Predicting radiation characteristics of a large array can be an 
easy task if the mutual coupling between the elements can be 
ignored. Pattern multiplication method of the array factor, and 
the single element radiation pattern of an isolated element or in 
the array environment, could be sufficient to predict the array 
radiation characteristics. However, determining the effective 
input impedances at the antenna terminals (input impedance 
while all the elements are active), is not a simple task and it is 
time-consuming if one needs to measure them. Since the demand 
for phased arrays, which provide electronic scanning capabilities 
has increased for satellite, military and vehicular 
communications increased, needs for accurate analysis and 
measurements becomes necessary. However, accurate design 
requires a full-wave analysis of the array. With the 
advancements in the computation methods, array analysis and 
design of large array have now become possible. However, such 
analyses can be highly computer-intensive and costly, since they 
require the use of expensive commercial codes [1, 2]. In fact, for 
very large arrays, such analysis may not even be feasible, unless we 
resort to parallel computing or use special numerical techniques. 


Although numerous design techniques can be found in the 
literature [3]-[22], these techniques are limited in their 
applicability to certain types of arrays with specific types of 
elements. Furthermore, it may be necessary to tailor the 
algorithm to the type of elements and, hence, modify it if 
element-type is changed. In some cases the formulation of the 
problem changes as the element type changes. The need for an 
efficient method to analyze large array is growing. 


In phased arrays, the mutual coupling has a pronounced 
effect on the performance as the scan range is increased. 
Mutual coupling changes the performance of a single element 
in the sense that the coupling between different elements is 
likely to affect the resonance frequency and the radiated 
power of each element. To account for this, there has been a 
considerable interest in developing techniques for including 
the mutual coupling effects into the design techniques for 
phased arrays; however, such techniques are both complex 
and computer intensive. To reduce the design costs and risks, 
and to improve the performance of the arrays simulations, the 
design method should be computationally efficient, include 
effects of mutual coupling, and should be able to determine 
the array performance for the entire scan range. Furthermore, 
it should be able to predict the matching performance within 
the scan range and the effect of failing elements in the array. 
Therefore, in order to design these arrays and their feeding 
network to accommodate the wide-angle scanning capabilities 
and to account for their effect on the matching of the array 
elements, a computationally efficient method needs to be 
developed. Usually, for a very large finite array, the analysis 
of single element in an infinite array environment is 
performed numerically, and the analysis of the large arrays is 
then extrapolated from this analysis. However, the elements 
near the edge perform differently from the rest of the array 
elements, and the analysis may have to be repeated, for each 
scan direction. Simulations based on the generally employed 
infinite-array approach and simulations based on the periodic 
version of the finite-element method do not satisfy these 
criteria. Simulations using the infinite array environment do 
not describe the edge effects, while simulations of large finite 
arrays are computationally too expensive. Furthermore, 
neither of these simulations provide direct insight into the 
physics relevant to the design. To overcome these 
disadvantages, simple simulation tools for large arrays are 
desired. 
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In this paper, the technique described in [23, 24] is used to 
develop a procedure for the design of planar antenna arrays. The 
design expressions are used in a rather general form to analyze 
arrays of similar elements. This technique has been presented in 
[25] and implemented in a near-field focusing array [26]. Here, 
we propose a technique which leverages the advantages of 
numerical techniques and the availability of efficient 
commercial codes, such as WIPL-D [27], to analyze a small 
array, and use the computed parameters to predict and design a 
larger array, which has a similar lattice. 


II. ANALYSIS OF PLANAR ARRAY 


For an array of N elements, the antenna element is 
considered as a load in series with internal resistance of the 
source as shown in Fig. 1. The mutual coupling between the 
elements changes the elements terminal voltages as shown in 
Fig. 2. It is important to derive a relationship between the 
effective voltage across the input terminals of the antenna and 
the applied voltage (source voltage). If J; is the current 
supplied by the source voltage V; to the i-th radiating element, 
the effective voltage V4; at the input terminals of that element 
is determined through the relation [23] 
T; = Vj -Vai) l Ri (1) 


where R; is the internal resistance of the source and i=1,2...N.. 
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Fig. 1 Equivalent circuit of an antenna array 
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Fig. 2 Equivalent circuit for the effective voltages of the antenna array. 


Fig. 2 indicates that J; is linearly related to the effective 
voltages across input terminals. 


liE YV y (2) 


where Y; is the self-admittance of the i-th element and Y; is 
the mutual admittance seen by the feed between the i-th and j- 
th radiating elements. From (1) and (2), the effective voltage 
can be expressed as 


ra VG Ri SV yV: (3) 


where 6; = 1 for i = j and zero for i #j. This equation can be 
expressed in a matrix form as 


2 


[Val=[Ai;] LV ] 


where 


-1 
[Aj ]5LY:; Ritó] (5) 
The effective voltage across the i-th element can be derived as 


Va ZE AV; (6) 


The effective voltage at the input terminals of each element 
can be obtained from the known source distribution. If the 
mutual coupling effects are not included, then the applied 
voltage V; appears directly across the input terminals of the 
patch. 


(4) 


Effective Input Impedance: The effect of mutual coupling 
on the input impedance of each element in an array can be 
expressed as 


Zi Ve ‘Ti. 


(7) 
which in turn, can be rewritten as 
Z= ae Ai; V ;/ ee Yir Ag V j (8) 
Also, using (1) to (4) Zi can be expressed as 
Zi =R; E Aij Vif Etal Si- Ag )V y (9) 


for R;+0. Ifthe mutual coupling effects are not included, then 
the input impedance of the element is equal to that of an 
isolated element. It should be mentioned that the mutual 
impedance matrix is independent of the excitation sources that 
are used to control the scanning and the side lobe level. 
However, the effective impedances strongly depend on the 
mutual impedance matrix and the excitation of the array. This 
simple procedure could also be used to study the effect of 
shorting or open circuiting any antenna port within the array 
by changing the value of the internal impedance of the source. 


Radiation Patterns of a Planar Array: Consider an array of 
N identical elements in xy-plane. The total field pattern of 
such an array is the product of the element pattern and the 
array factor. If (p;, @;) are the cylindrical coordinates of the i- 
th element and £; is the effective excitation voltage, then the 
array factor is expressed as 


FOD EG erred (10) 


Then, the expressions for total pattern of the array of 
identical elements are given by 


Eb Eo 
= F(0,¢) 


where Eo and E, are the element patterns, within the array 
environment. For a phased array, the main beam direction is 
controlled by the phase of the excitation sources. The phase 
shift, which should be introduced in the excitation voltage of 
the i-th element to obtain a main beam in the direction of (00, 
Qo), is given by 


W,=- jk P;8in 8, cos (¢,- ¢,) (12) 


where ¢; is the phase of the excitation voltage V;. However, 
the mutual coupling affects the effective voltage magnitude 
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and phase across the antenna terminals. Therefore, it is 
possible that the main beam direction and the desired side lobe 
levels be altered. 


Power Gain of an Array Antenna: With the knowledge of 
effective input impedance and the excitation current, the 
power delivered to an element can be obtained, where the real 
part of the power delivered to the i-th element, which 
contributes to the radiation is given by 


e 2 
P =R; |1| (13) 
where Rin’ is the real part of the effective input impedance Z£. 
The total input power that contributes to the radiation from the 
array is obtained as 
N N 


R= 2) KIL) 


i=l i= 


= 


; (14) 
If the radiated power is assumed to be the power radiated 
by an isotropic source, the power radiated per unit area by the 
isotropic source is Pr= P7/(4nr?). The total radiated power 
density available from the array at a distance r is given as 


“LEP +1 Eye 
array 1207 (15) 


The gain of the array with respect to an isotropic radiator can 
be expressed as 


| £5 È +E? 


GO, ) = 4ar 
= Ep 


(16) 
Construction of the Mutual Admittance Matrix for a 
Large Array: 


For simplicity, we will illustrate the procedure, proposed 
here in, by using a 4-element linear array, with N, = 4, to 
construct a linear array of 8 elements, as shown in Fig. 3. The 
mutual admittance matrix of the small array is of order (4x4) 
can be presented graphically as shown in Fig. 4a. We simply 
map the small matrix into the large matrix for the 8-element 
array, which is an 8x8 matrix, as shown in Fig. 4b. Note that 
the mutual coupling between elements with distances greater 
than 3d are ignored and are forced to be zero. The mutual 
couplings between the elements of the same relative distances 
and position in the large matrix are retained the same as those 
in the small matrix. It is also advisable to keep track of the 
relative positions in the small and large array. 


z 
ALLELDI 


Fig. 3. Linear array of 4 elements that is used to construct a linear array of 
8-elements. 
































1,1 1,2 1,3 1,4 

2,1 22 2,3 2,4 

3,1 3,2 3,3 3,4 

4,1 4,2 4,3 4,4 
(a) Yj; (for the 4-element array) 
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(b) Yj; (for the 8-element array) 


Fig. 4 Mutual admittance matrix arrangement (a) 4-elemet array and (b) 
constructed for 8-element array. 

For the two-dimensional array, similar logic is used to 
construct the mutual impedance matrix of the larger array. The 
large array is represented by Mz elements in the x-direction 
and Nz in the y-direction. The small array has M, and N, 
elements in the x-and y-directions, respectively. Notice that 
Mı > M; and Nz > Ns. If we denote the admittance matrix 
elements of the large array by the superscript L and those for 
the small array by the superscript S, then we have 

S roa 
yt- Y; i, j =1,2,.., M N, 

ke T . 
0 otherwise 
(17) 


where 
bist, =1,2,..., M;N; 
k =M,;(j; Slag 
i=M (j -D +i, 


l=N,0,-)+J,' 
IAN GAUSS! 
and 


i, =i =li =i &l j Hla. Fe | 


HI. NUMERICAL RESULTS 


To test the procedure described in the last section, an array 
of rectangular microstrp patches excited by a hook-shaped 
probe is considered, as shown in Fig. 5. The patch has the 
dimensions: length 25 mm; width 30 mm; and the height 
above the ground plane 7.5 mm. The vertical part of the 
probe that is above the ground plane and just below the patch 
edge is 6 mm. It is bent horizontally and extended under the 
patch by 11 mm and then bent vertically toward the ground 
plane by 2 mm. The ground plane is assumed to be infinite. 
The input impedance and radiation patterns of this patch are 
shown in Figs. 6 and 7, respectively. The radiation patterns 
are obtained at 4.5 GHz. The results are obtained by using 
the commercial code WIPL-D [25], which is based on the 
surface integral equations and the Method of Moment 
(MoM). 
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Also, this array is analyzed using the WIPL-D code to obtain 
the effective impedance of the array. Fig. 8 shows a 
comparison between the effective impedances using the 
MoM and the one obtained with the present method. Notice 
that the array was numbered row by row. It can be observed 
that the edge elements are represented by the effective 
impedance spikes in both the real and imaginary parts. The 
radiation patterns of this array, shown in Fig. 9, using array 
pattern multiplication (APM) compared with the MoM, 
which include the mutual coupling effects. It should be 
mentioned that the APM method used the radiation pattern 
of an isolated element. If the radiation pattern of the element 
in the array environment is used, an improvement can be 
obtained particularly at the lower elevation angles. The array 
gain is computed and found to be 25.86 dB using the APM, 





Fig. 5 Geometry of a rectangular patch excited by a Hook-shaped probe. 
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Fig. 6 Input Impedance of the rectangular patch 
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Fig. 8 Effective impedance of 8x8-element array using the 4x4-element 
array. The dotted curves are the actual impedances from the full wave 
analysis (WIPL-D). 




















(b) E-plane 
Fig. 7 Radiation patterns of the patch at 4.5 GHz. 
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A 4x4-element array using the patch described above, 
with center-to-center separation distance of 38 mm between 
the elements in the x-direction and the y-direction, is 
analyzed using the WIPL-D code. The code provides the -40 
mutual admittance matrix, whose size is on the order 
16. This matrix is used to fill the mutual admittance matrix 
of an 8x8-element array. The effective impedance of the 


8x8-element array is computed using the proposed Fig. 9 Radiation patterns of the 8x8-element array. The solid and dashed 
procedure. lines are for the MoM code. 
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Next, a larger array of 16x16-element array is constructed 
by using the 4x4-element array. The effective impedances 
are shown in Fig. 10. The radiation patterns are plotted in 
Fig. 11. The computed gain is found to be 31.6 dB when 
using the APM vs. 31.5 dB using the MoM. 
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Fig. 10 Effective impedance of 16x16-element array using the 4x4- 
element array. 
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Fig. 11 Radiation patterns of the 16x16-element array using 4x4- 
element array. 

When this array is constructed from an 8x8-element array, 
the effective input impedance is given in Fig. 12. It can be 
observed that the impedance prediction becomes more 
accurate than when the 4x4-element array is used. Fig. 13 
shows the radiation patterns. Comparing the radiation 
patterns in Fig. 11 with those in Fig. 13, we observe better 
prediction of the far away side lobe levels. 
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Fig. 12 Effective impedance of 16x16-element array using the 8x8- 
element array. 
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Fig. 13 Radiation patterns of the 16x16-element array using 8X8- 
element array. 

Finally, we move to a 32x32-element array and predict its 
performance by using a 4x4-element array; an 8x8-element 
array; and a 16x16-element array. With the current resources 
we were not able to predict one MoM results for the 32x32- 
element array. The effective impedance for the 32x32- 
element array is shown in Fig. 14. From the figure one can 
see that the 8 x 8-element arrays provide results that are 
similar to those of the 16x16-element array. Also observe 
that the 4x4-element array provides a reasonably good 
prediction for the 32x32-element array and that the results 
predicted by the 16x16-element array are virtually the same 
as those for the 32x32-element array (not obtained from full- 


wave analysis). The 32x32-element array achieves a gain of 
37.46 dBi. 
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Fig. 14 Effective impedance of 32x32-element array using 4x4 (black), 8x8 
(yellow), and 16x16-element array (blue). 

Using the same mutual impedance values, we can use 
different excitations to derive the effective voltages and 
currents and compute the effective impedances. That is to 
say using the same A matrix in (4), we can change the vector 
|V> to obtain different array characteristics. Such 
information can be used to predict the array performance 
under different variations of scanning conditions and to 
know the range of the effective input impedances of the array 
elements. To illustrate that, let us consider the case of a 
linear array of 32 elements that is obtained from a 4 element 
array, whose center-to-center separation is 38 mm. Figs.15- 
16 show the effective impedance of the array for the H-plane 
(hook probe plane is normal to the array line) and E-plane 
(hook probe plane is parallel to the array) arrangements, 
respectively. The antenna gains were predicted to be 23.7dBi 
and 23dBi, respectively. When we scan the array off the 
broadside by 45°, the input effective impedances change as 
shown in Figs 17 and 18. The radiation patterns predicted by 
the AFM are in excellent agreement with the actual antenna 
predicted by the MoM. The gain of the antenna is predicted 
to be 23.0 and 22.5 dBi, respectively. 


125 





100 $ 





O-O-9-0-0-0-Ox.4 


TO 





present 
o actual 





Ci 

pe 
50 

BB 

s 
25 + 


OF x 





O-0-9-0-0-0-0-0- 








-25 1 1 1 
8 16 24 32 


element number 
Fig. 15 Effective input impedance of an H-plane 32-element linear array 
with uniform distribution obtained from 8-element array and compared 
with the actual array. 
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Fig. 16 Effective input impedance of an E-plane linear array of 32 
elements with uniform distribution obtained from 8-element array and 
compared with the actual array. 
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Fig. 17 Effective input impedance of an H-plane linear array of 32 
elements with uniform distribution and scanning at 45° and compared with 
the actual array results. 

When we use the Taylor distribution to reduce the 
sidelobe level to -30 dB, the effective input impedances are 
obtained as shown in Fig. 19 for the 32-element E-plane 
array. It can be seen that the impedance level is not that much 
different from those plotted in Fig.17. The radiation pattern 
is shown in Fig. 20. When this array, designed to meet the - 
30 dB sidelobe level requirement is scanned at 45° the 
effective input impedance is shown in Fig. 21 and the 
radiation patterns are shown in Fig. 22. A significant change 
in the effective impedance behaviour is observed and the 
radiation patterns show a grating lobe at 90°. The mutual 
coupling pushes the sidelobe above the level of -30dB. It 
should be mentioned that the radiation patterns for the 
uniform power distribution case does not show the grating 
lobe effect since it is below the regular side lobe level. 
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Fig. 18 Effective input impedance of an E-plane linear array of 32 
elements with uniform distribution and scanning at 45° and compared with 
the actual array results. 
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Fig. 19 Effective input impedance of an E-plane linear array of 32 
elements with Taylor distribution of -30dB sidelobe level compared with 
the actual array results. 
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Fig. 20 Radiation patterns of the array in 19. 








500 
400 L pore 
ooo? 
300 | gaeth 
Ci pA eee 
s 2-0-0-0 present 
P< ooo L ° o actual 
R oOo 
-F OPCS 006 4 g 
PA 100 | i 
Ooo O 
o L 
-100 











8 16 24 32 
element number 


Fig. 21 Effective input impedance of an E-plane linear array of 32 elements 
with Taylor distribution and -30 DB sidelobe level and scanning at 450 
compared with the actual array results. 
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Fig. 22 Radiation patterns of the array in Fig. 21. 


The method described above can also be used to study the 
effect of faulty elements, such as shorted or open circuited 
elements, or combinations thereof. However, the matrix A 
has to be constructed using (5) based on the mutual 
admittance Y of the array and the new internal impedances 
of the sources. 


IV. CONCLUSIONS 


A simple and efficient procedure has been proposed to 
predict the performance of large finite arrays from the 
measurements results of a smaller array. It is necessary to 
ensure that both arrays have the same lattice arrangements and 
spacing between the elements. The analyses have been 
performed numerically; however, experimental study can be 
used with the small array and the large finite array 
performance can be predicted based on the measurements of 
the mutual impedances of the small array. The procedure has 
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been tested by using an array of air-filled microstrip patch 
antennas. The results were verified to provide confidence on 
the proposed procedure. A future analysis will be presented 
based on the S-parameters of the array and its frequency 
response, different element types will be investigated and the 
results will be reported in a future publication. 
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